The significance of the fragmentation region in ultrarelativistic heavy ion collisions 
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We present measurements of the pseudorapidity distri- 
bution of primary charged particles produced in Au + Au 
collisions at three energies, y?s NN = 19.6, 130, and 200 GeV, 
for a range of collision centralities. The centrality dependence 
is shown to be non-trivial: the distribution narrows for more 
central collisions and excess particles are produced at high 
pseudorapidity in peripheral collisions. For a given centrality, 
however, the distributions are found to scale with energy 
according to the "limiting fragmentation" hypothesis. The 
universal fragmentation region described by this scaling grows 
in pseudorapidity with increasing collision energy, extending 
well away from the beam rapidity and covering more than half 
of the pseudorapidity range over which particles are produced. 
This approach to a universal limiting curve appears to be 
a dominant feature of the pseudorapidity distribution and 
therefore of the total particle production in these collisions. 



The strong interaction, described by Quantum Chro- 
modynamics (QCD), may be studied under conditions 
of high parton density and high energy density, using 
ultrarelativistic heavy ion collisions. The high den- 
sity regime of QCD is sensitive to nonlinear dynamics 
and nonperturbative effects, including parton saturation, 
the onset of color deconfmement and chiral symmetry 
restoration. More specifically, the pseudorapidity density 
of charged particles dN c h/dr], where i] = — lntan((9/2), is 
related to the entropy density at freezeout. It has been 
demonstrated that the growth with energy of dN c h/dr] at 
mid-rapidity is modest compared to the original expecta- 
tions [1,2], and provides a strong constraint on the initial 
state parton density and further particle production dur- 
ing the subsequent evolution of the system. This Letter 



focuses on particle production away from mid-rapidity, 
which constrains the collision dynamics more completely. 

We have measured the pseudorapidity distribution of 
charged particles, dN c h/dr], over a broad range of r) for 
Au + Au collisions at a variety of collision centralities 
(impact parameters). These measurements were made 
for three energies, t/s^ — 19.6, 130, and 200 GeV, 
covering a span of an order of magnitude in the same 
detector, allowing for a reliable systematic study of par- 
ticle production with energy in these collisions. The data 
were taken using the PHOBOS apparatus [3,4] during the 
year 2000 and year 2001 runs of the Relativistic Heavy 
Ion Collider at Brookhaven National Laboratory. The 
apparatus used in this analysis comprises a set of silicon 
detectors with nearly complete coverage for < 5.4, 
which are used for detecting the charged particles, and 
plastic scintillator counters, covering 3 < \rj\ < 4.5, used 
for triggering. Event selection cuts were made on the 
plastic scintillators or silicon detectors and also on tim- 
ing information from the forward hadronic calorimeters, 
which measure the energy deposited by spectator neu- 
trons. 

The pseudorapidity densities dN c h/drj were corrected 
for particles which were absorbed or produced in dead 
material (primarily the beryllium beampipe or the mag- 
net steel) and for feed-down products from weak decays 
of neutral strange particles. More details of the analysis 
procedures leading to the charged particle pseudorapidity 
density can be found in Ref. [5]. Two improvements in 
the handling of common mode noise in the silicon detec- 
tors have been implemented for this analysis, leading to 
small changes in the results and to a slight reduction in 
the systematic errors. First, the event-by-event common- 
mode noise in the ring detectors (3.0 < |7y| < 5.4) was 
found to grow with pad size, and the correction scheme 



was modified to include this effect. This refinement 
was already in place in Ref. [2]. Second, for very high 
occupancies, the common mode noise correction in the 
octagon detector (\rj\ < 3.2) becomes slightly inaccurate. 
A comparison of the data from the octagon detector with 
those from the more highly segmented vertex detector 
was used to determine a correction factor for this effect 
as a function of 77, centrality, and beam energy. This 
correction was only required near mid-rapidity (| r\ \ < 1.5) 
for the central data and it was less than 4% everywhere 
even in the 0-6% centrality bin of the 200 GeV data. For 
the 6-15% bin of the 200 GeV data and the 0-6% bin 
of the 130 GeV data, this correction was less than 1% 
everywhere. For the more peripheral high energy bins 
and for the 19.6 GeV data, no correction was required. 

The centrality of the collision is characterized by the 
average number of nucleon participants (N part ). For the 
130 and 200 GeV data sets, this was estimated from 
the data using two sets of 16 paddle counters covering 
3 < < 4.5 forward and backward of the interaction 
point. The truncated mean of the signals in each set 
of detectors is proportional to the total charged mul- 
tiplicity in this region of pseudorapidity. In order to 
extract (N part ) for a given fraction of the cross-section, 
we rely on the fact that the multiplicity in the paddles 
increases monotonically with centrality, but we do not 
assume that the number of participants is proportional 
to the multiplicity. The monotonic relationship between 
the truncated paddle mean and N part was verified us- 
ing the neutral spectator energy measured in the for- 
ward hadronic calorimeters. For peripheral events, the 
dominant uncertainty in (N part ) is given by the trigger 
efficiency uncertainty. Using several methods, based on 
the HIJING model [6] , we have estimated our minimum- 
bias trigger efficiency for events with a vertex z vtx near 
the nominal interaction point to be e = 97 ± 3%. The 
final systematic errors on (N part ), ranging from 3-6%, 
are tabulated in Ref. [7] where this method was outlined 
in more detail. 

At the lowest RHIC energy, ^fs^ = 19.6 GeV, the 
much lower beam rapidity (jjbeam ~ 3) precludes us- 
ing the same method of analyzing the trigger counters 
(3 < \r)\ < 4.5). Instead, we construct a different quan- 
tity, "EOCT" , which is approximately proportional to 
the multiplicity: the path-length corrected sum of the 
energy deposited in the octagon (silicon) detector(|?7| < 
3.2). In the first step, both HIJING simulations and 
simple Glauber model [8,9] Monte Carlo simulations are 
used to estimate the fraction of the total cross-section (e) 
in the triggered sample, as well as to estimate the system- 
atic error. For consistency with the higher energy results, 
the low N par t tail is assumed to be distributed as in 
HIJING, leading to e = 52 ± 4%. Once e is determined, 
cuts are made in EOCT in a similar way as with the 
paddle signal to extract (N part ) for a chosen fraction of 
the total cross-section [7] . As with the paddle signal, we 
only assume that EOCT is monotonic with (N part ), not 
that it is proportional to it. The (N part ) values and the 



estimated systematic uncertainty for various centrality 
bins in the 19.6 GeV data are given in Table I. 
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FIG. 1. The charged particle pseudorapidity distribution, 
dN c h/dri, measured for Au + Au at ^/s NN = 200, 130, and 
19.6 GeV for the specified centrality bins. These bins range 
from 0-6% central to 45-55% in the case of the higher energy 
data and 0-6% to 35-45% for the 19.6 GeV data. The 
statistical errors are negligible. The typical systematic errors 
(90% C.L.) are shown as bands for selected centrality bins. 

Figure 1 shows the charged particle pseudorapidity 
distributions (dN c h/drj) measured at ^fs NN — 200, 130, 
and 19.6 GeV for different centrality bins for —5.4 < 
r\ < 5.4. The statistical errors are comparable to the size 



of the data points. The 90% confidence-level systematic 
errors are shown as a gray band for selected, but typical, 
centrality bins. 

Due to the large coverage in 77, dN ch /dri is measured 
over almost the full range, except for a small missing 
fraction at very high \rj\. Using data at lower energy 
(see Figure 2), this fraction is estimated to be less than 
2% for the 130 and 200 GeV central (0-6%) results. 
Table II lists the integrated charged multiplicity for three 
energies in different ranges of 77. The fiducial range 
1 77 1 < 4.7 is quoted for easier comparison to data from 
the BRAHMS experiment [10,11]. After correcting for 
the slight differences in centrality binning, the BRAHMS 
data for 130 and 200 GeV are each ~ 6% below the results 
quoted here, consistent within the quoted errors of both 
experiments. The remaining columns show the result in 
the fiducial range |ry| < 5.4 of PHOBOS, and the total 
inferred N c h- 

To separate the trivial kinematic broadening of the 
dN c h/drj distribution from the more interesting dynam- 
ics, the data for Au + Au collisions at different energies 
can be viewed in the rest frame of one of the colliding 
nuclei. Such an approach led to the ansatz of "lim- 
iting fragmentation" [12], which successfully predicted 
the energy dependence of particle production away from 
mid-rapidity in hadron collisions, including pp [13]. This 
ansatz states that, at high enough collision energy, both 
d 2 N/dy'dp T and the mix of particle species (and there- 
fore also dN/dr)') reach a limiting value and become 
independent of energy in a region around y' *~ 0, where 
y' = y - y b earn and rapidity y = tanlT 1 f3 z . 

Figure 2a shows the scaled, shifted pseudorapidity 
distributions dN c h / drj' / (N part / 2) for central Au + Au 
collisions which span a factor of ten in collision energy 
(V s «jv)- The distributions are scaled by (N part /2) to 
remove the effect of the different number of nucleons 
participating in collisions with different centralities. The 
results are folded about mid-rapidity (positive and nega- 
tive 77 bins are averaged) . The distributions are observed 
to be independent of collision energy over a substantial 77' 
range. This is consistent with and extends a similar ob- 
servation made by BRAHMS [11] over a more restricted 
77' range. Both the 19.6 and 130 GeV data reach 85-90% 
of their maximum value before deviating significantly 
(more than 5%) from the common limiting curve. These 
data demonstrate that limiting fragmentation applies 
well in the Au + Au system and that the "fragmentation 
region" is rather broad, covering more than half of the 
available range of 77' over which particles are produced. In 
particular, the fragmentation region grows significantly 
between 19.6 GeV and 130 GeV, extending more than 
two units away from the beam rapidity. A similar effect 
was observed in pp collisions, over a range of a/s from 
53 to 900 GeV [13]. In both cases, particle production 
appears to approach a fixed limiting curve which extends 
far from the original beam rapidity, indicating that this 
universal curve is an important feature of the overall 
interaction and not simply a nuclear breakup effect. 



This result is in sharp contrast to the boost-invariance 
scenario [14] which predicts a fixed fragmentation region 
and a broad central rapidity plateau that grows in extent 
with increasing energy. 
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FIG. 2. Au + Au data for v /s^7 = 19.6, 130, and 200 GeV, 
plotted as dN c h/dr)' per participant pair, where r/ = r\ — yteam 
for a) 0-6% central, b) 35-40% central. Systematic errors 
(90% C.L.) are shown for selected, typical, points. 

Figure 2b shows the scaled pseudorapidity distribu- 
tions for a set of non-central collisions, which also exhibit 
limiting fragmentation over a broad range of 77'. These 
results also show that the extent of the fragmentation 
region is nearly independent of centrality 
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FIG. 3. The distribution dN c h/drj' per participant pair for 
central (0-6%) and non-central (35-40%) Au + Au collisions 
for a) = 200 GeV, b) ^/s^ = 19.6 GeV. Systematic 

errors are not shown. 
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FIG. 4. The ratio of dN c h/drj' per participant pair between 
non-central (35-40%) and central (0-6%) data plotted for 
v /s^T = 200 GeV, 130 GeV, and 19.6 GeV. The errors 
represent a 90% C.L. systematic error on the ratio. 

Figure 3 shows the centrality dependence of the 
dN/drf /(Np ar t/2) distribution at the two extreme ener- 
gies: 19.6 and 200 GeV. These data demonstrate that the 
result in the fragmentation region changes significantly 
with centrality. Figure 4 shows the ratio of non-central to 
central data with (90% C.L.) systematic errors included. 
The error in the ratio involves a partial cancellation of the 
systematic errors in the individual measurements. For 
rj > —1.5, the scaled pseudorapidity density actually 
grows in the peripheral data with respect to the more 
central data. This effect has already been observed for 
Au + Au collisions at y/s^ — 130 GeV [5] and for 
Pb + Pb collisions at 17 GeV [15], and is confirmed here 
for Au + Au at 19.6 and 200 GeV. This contradicts the 
suggestion, put forward in references [10,11], that the 
limiting curve for particle production in the fragmen- 
tation region is independent of centrality as well as of 
energy. It should be noted that the hypothesis of limiting 
fragmentation does not imply that the limiting curve is 
independent of centrality, just that dN / dij' / (N part / 2) is 
energy-independent for a fixed centrality, once at suffi- 
ciently high energy. 

The strong centrality dependence seen in Figure 3 
has two features: an excess of particles at high rj in 
peripheral Au + Au events and a narrowing of the overall 
pseudorapidity distribution, perhaps caused by a shift 
of particles from high rj toward mid-rapidity for more 
central collisions. The excess of particles at high \rf\ can 
be most easily seen in Figure 1, where, for the 19.6 GeV 
data, the absolute yield of charged particles actually 
grows for non-central collisions. It should be noted that 



the very peripheral limit of a Au + Au collision is not an 
N + N collision, but rather an N + N collision with two 
large excited nuclear remnants, one of which is at y' — 0. 
The increased particle production near y ± yteam = 
due to the nuclear remnants has been studied in p + A 
collisions [16] and for lower energy Pb+Pb collisions [15]. 
The narrowing of the overall distribution may be due 
to dynamical effects, such as baryon stopping [17], or 
kinematic effects, such as a shift in rj' (for fixed y') due 
to the particle mix (p/7T ratio) changing with centrality. 
This narrowing of the distribution can be characterized as 
an increase at mid-rapidity [7,18] with an approximately 
compensating reduction at high rj. 

In summary, we have performed a comprehensive ex- 
amination of the pseudorapidity distributions of charged 
particles produced in Au+Au collisions at RHIC energies 
from y/s^" = 19.6 to 200 GeV, including an estimate of 
the full charged particle multiplicity at three energies. 
For central collisions at the highest energy, we find that 
a total of more than 5000 charged particles are produced. 
These results span eleven units of pseudorapidity, a factor 
of ten in energy, and a factor of five in (N part ) — all 
measured in a single detector. The data show a number 
of interesting features. First, limiting fragmentation 
(energy independence of dN / drj 1 ) is valid over a large 
range of rj' . Second, the scaled dN/drj shape is not 
independent of centrality at high 77. The rj distribution is 
broader in peripheral collisions than in central collisions. 
Third, as in pp collisions, the fragmentation region in 
Au + Au collisions grows in pseudorapidity extent with 
beam energy, becoming a dominant feature of the pseu- 
dorapidity distributions at high energy. 
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TABLE I. Estimated number of nucleon participants in 
19.6 GeV Au + Au collisions according to the centrality bin 
(quoted as a fraction of the inelastic cross-section). The 
systematic errors are shown. Note: the 40-45% centrality 
bin has large systematic errors on {N pa rt}, so N par t scaling is 
not used for that bin. 



Centrality {N par t} 




Centrality 


{Np ar t} 


0-3% 


351 ± 14 




20-25%: 


178 ± 12 


3-6% 


322 ± 10 




25-30%: 


150 ± 13 


6-10% 


287 ± 10 




30-35%: 


125 ±13 


10-15% 


247 ±11 




35-40%: 


103 ± 14 


15-20% 


210 ±12 








0-6% 


337 ± 12 




15-25%: 


194 ± 12 


6-15% 


265 ±11 




25-35%: 


138 ± 13 


TABLE II. Total charg 


;cd 


multiplicity in 


three fiducial 


ranges of 


r\ for central (0-6%) 


collisions. 




\/ S NN 


N ch (\ V \ < 4.7) 


N ch (\rj\ < 5.4) 


N ch (total) 


19.6 GeV 


1670 ± 100 


1680 ± 100 


1680 ± 100 


130 GeV 


4020 ± 200 


4100 ± 210 


4170 ± 210 


200 GeV 


4810 ± 240 


4960 ± 250 


5060 ± 250 



